A novel solution to optimal power flow problems using composite differential evolution integrating effective constrained handling techniques

Optimal power flow is a complex and highly non-linear problem in which steady-state parameters are needed to find a network’s efficient and economical operation. In addition, the difficulty of the Optimal power flow problem becomes enlarged when new constraints are added, and it is also a challenging task for the power system operator to solve the constrained Optimal power flow problems efficiently. Therefore, this paper presents a constrained composite differential evolution optimization algorithm to search for the optimum solution to Optimal power flow problems. In the last few decades, numerous evolutionary algorithm implementations have emerged due to their superiority in solving Optimal power flow problems while considering various objectives such as cost, emission, power loss, etc. evolutionary algorithms effectively explore the solution space unconstrainedly, often employing the static penalty function approach to address the constraints and find solutions for constrained Optimal power flow problems. It is a drawback that combining evolutionary algorithms and the penalty function approach requires several penalty parameters to search the feasible space and discard the infeasible solutions. The proposed a constrained composite differential evolution algorithm combines two effective constraint handling techniques, such as feasibility rule and ɛ constraint methods, to search in the feasible space. The proposed approaches are recognized on IEEE 30, 57, and 118-bus standard test systems considering 16 study events of single and multi-objective optimization functions. Ultimately, simulation results are examined and compared with the many recently published techniques of Optimal power flow solutions owing to show the usefulness and performance of the proposed a constrained composite differential evolution algorithm.

• Two representative constraint techniques, such as feasibility rule (FR) and epsilon constraint method (ECM), and their combinations are employed with the current state-of-the-art unconstrained CoDE search algorithm to solve the OPF problem.• Sixteen events of highly complex non-linear objective functions are formulated to solve single and multi- objective OPF problems and show the superiority and performance of the proposed algorithm.• Simulation results of all the algorithms C2oDE-FR, C2oDE-ECM, C2oDE-FR-ECM, and C2oDE-ECM-FR are statistically compared.• Small to large-scale power system networks such as IEEE 30, 57, and 118-bus networks are adopted to test the proposed Algorithm.
The remaining division of this article is planned as Sect. 2 contains mathematical modeling of OPF and constraint handling techniques, and Sect. 3 describes the objective function and study events.The proposed optimization algorithm is defined in Sect.4, simulation results and comparisons are discussed in Sect.5, and concluding remarks are produced in Sect.6.

Mathematical modelling of OPF problem
Generally, OPF is a complex and non-linear problem, and its main objective is to optimize single and multiobjective functions subject to satisfy the set of equality and inequality constraints.Mathematically, the OPF problem is described as follows: Minimize

Constraints
The solution to the OPF problem must achieve both equality (active and reactive power balance) and inequality (operating limits of power system components) constraints.Figure 2 shows the equality and inequality constraints examined in the present study.
where, P D i and Q D i are the active and reactive demand at bus i, G ij and B ij are shunt conductance and susceptance between bus i and j respectively.δij is the voltage angle difference between bus i and j and shows NB the number of buses.
At the time of the optimization process, the proposed algorithm chooses the values of each variable between the min and max limit.

Proposed constraint handling techniques
Usually, all the real word problems are constraint type defined in Eq. (1), in which the equality constraints h j ( − → x , − → u ) given in Eqs. ( 2) and ( 3) are automatically satisfied when the solution of power flow is converged.However, special attention is needed to inequality constraints g j − → x , − → u given in Eqs. ( 4) to (10).Generally, the j th inequality constraint violation G j ( − → x ) is given as: However, the overall degree of constraint violation G( − → x ) can be calculated by the sum of all the inequality constraint violations and given as: Constrained optimization problems mean to search in the feasible region, and EAs are population-based stochastic search methods in which an infeasible solution is complicated to discard.Therefore, proper CHTs are used together with EAs to enhance the overall performance of an algorithm.This work proposes two CHTs: feasibility rule (FR) and ε constrained method (ECM).FR is given in 22 and suggests three rules to compare any two solutions described as follows: 1.Both solutions are feasible; select the one with a better objective function value.2. Both solutions are infeasible; choose the one with a lower value of constraint violation.3.One is feasible, and the other is infeasible; always select a feasible one.
The second proposed CHT is ECM has given in 23,24 , in which two solutions − → x i and − → x j are compared as follows: In (13), the parameter ε decays as increasing the iteration number and is given as where the parameter ε 0 is the primary threshold, initially it is equal to max(G − → x ) , T is the maximum generation, t is the current generation, constant parameter λ = 6 recommended in 25 and p controls the degree of convergence of objective function.

Objective functions and study events
To highlight the superiority and effectiveness of the proposed C2oDE algorithm by considering the various CHTs, 16 events comprised of single and multi-objective functions are evaluated and implemented on IEEE 30, 57, and 118-bus standard IEEE networks.Bus 1 is considered the slack/reference bus in the event of 30 and 57-bus systems; however, in the 118-bus system, the 69th bus is the slack/reference bus.The role of the reference bus is (5) www.nature.com/scientificreports/ to achieve equality constraints given in Eqs. ( 2) and ( 3) during the load flow study.In subsequent sub-sections, the mathematical formulation of different events for the 30, 57, and 118-bus tests is described.

IEEE 30-bus system
The base MVA, bus, branch, and generator data of the IEEE 30-bus test network is taken from 26 , and a summary of the significant components of this system is arranged in Table 1.There are 10 events are formulated for the IEEE 30-bus network, in which the first six events comprised of minimizing single objective and the remaining four events are based on weighted sum multi-objective optimization.

Event 1: minimization of basic fuel cost
Almost in all the literature, minimization of fuel cost was considered, and the relationship between the generator output power (MW) and the fuel cost ($/h) is given by a quadratic curve described as: where, P G i is the generated output power of i th bus and a i , b i , c i The constant cost coefficients of that generator are given in 5,27 and classified as in Table 2.

Event 2. minimization of fuel cost multi-fuels
Thermal power generation may have multi-fuel resources, including coal, oil, and natural gas.Therefore, the relationship between fuel cost vs output power for such plants is given in the piecewise quadratic function shown in Fig. 3. Mathematically, the cost function of a multi-fuel i th generator is given as follows: where, P G i is the generator output power within the specified range of P min G ik , P max G ik and k is the fuel type.The total fuel cost of the objective function can be calculated using Eq.(18).
In this event, the multi-fuel cost is proposed for the two generators and range of output power (MW) with their coefficients given in 5 and shown in Table 3, whereas, the cost for the other four generators is identical as in the event 1.   (17).
The objective function of power system stability in this event is the maximum value of L j and is given as:

Event 4: minimization of emission
Many harmful gases such as SO x and NO x are emitted in tones per hour (t/h) into the atmosphere using conventional fuel's thermal power generation (MW).In the present event, the emission is considered the objective function of OPF and computed as: where, the values of the parameters α i , β i , γ i , ω i and µ i are given in Table 2.

Event 5: active power loss minimization
Mathematically active power loss (MW) can be given as: , where j = 1, 2, . . ., NL  where, G q(ij) is the conductance of branch q connected in between bus i and j and δ ij = δ i − δ j , is the voltage angle difference.

Event 6: minimization of basic fuel cost with valve-point loading
Valve-point loading wants to be measured for precise modeling and a more realistic cost of fuel vs generator output power (MW).Generation of power from multi-valve thermal engines shows variation in the fuel cost function, which is shown in the sinusoidal function.Such sinusoidal function is added to the fuel cost and resulting curve between output power (MW) vs fuel cost as shown in Fig. 4. Mathematically, generator fuel cost considering valve-point loading is given by 9 : where the constants d i and e i are the valve point loading parameters, and their values are given in Table 2.

Event 7: simultaneous optimization of basic fuel cost and active power loss
The weighted sum approach is used to convert multi-objective optimization functions into single-objective optimization and is denoted as: whereas, active power loss P loss can be computed bsing Eq. ( 23) and the parameter P is equal to 40 as suggested in 8 .

Event 8: simultaneous optimization of voltage deviation and fuel cost
According to power quality, the voltage deviation index is the most important aspect, and it is minimized by enhancing the voltage profile.The cumulative voltage deviation (VD) function at the PQ nodes is described as: The combined weighted sum of basic fuel cost and voltage deviation is given by: where the weight factor VD is assigned a value of 100 as in 9 and 8 .

Event 9: simultaneous optimization of voltage stability and fuel cost
Simultaneously, the minimization of basic fuel cost and maximization of voltage stability are converted into a single objective: whereas, the parameter L is called a weight factor equal to 100 suggested by 8 and L max is computed by Eq. ( 21).

Event 10: simultaneous optimization of cost, emission, losses, and vd
In this event, simultaneously, four objectives are considered to minimize, and the combined fitness function is given: where, E = 19, VD = 21 and P = 22 are the constant weights are considered the same as in 8 to balance among the objective functions.

IEEE 57-bus test system
To test the effectiveness of the C2oDE algorithm, the IEEE 57-bus system is considered.Four different events are considered to optimize with the C2oDE algorithm with two single objectives and the remaining two based on multi-objective, data given in Table 4.

Event 11: basic fuel cost minimization
In OPF, the basic objective is to minimize fuel cost, and mathematically, the function of fuel cost is the same as in Eq. ( 16).The coefficient of generator cost 26 and emission 5 are shown in Table 5.

Event 12: multi-objective optimization of fuel cost and vd
The weighted sum single objective optimization minimizes this event's basic fuel cost and VD.The fitness function in this study event is the same as in event 8 of IEEE 30-bus and mathematically is given by an Eq. ( 27).

Event 13: multi-objective optimization of voltage stability and fuel cost
The formulation of this event's weighted sum single objective function is the same as in event 9 of 30-bus.Also, lambda sub cap L is the same as in event 9.

Event 14: optimization of voltage deviation
In this event, the minimization of VD is considered the objective function of cumulative PQ buses and is calculated using Eq. ( 26).

IEEE 118-bus system
Furthermore, a large-scale 118-bus standard IEEE test network is considered to test the superiority of the proposed C2oDE algorithm.A couple of single objective events are considered for this system.Table 6 gives the bus, branch, generator, and other related data of the 118-bus network.
( The constant parameters of fuel cost are taken from 26 , and the formulation of the fuel cost function is similar to event 1 of 30-bus.

Event 16: active power loss minimization
In this event, the minimization of real power loss is considered the objective function and calculated using Eq. ( 23).

Proposed optimization algorithm
OPF is a constrained optimization problem, and how to solve constrained optimization problems has greater practical significance.Evolutionary algorithms (EAs) have involved noticeable attention in efficiently resolving practical constrained optimization problems in the past two decades.The constrained EAs have two main components: the search algorithm and the appropriate constrained handling method.Differential evolution (DE) is a popular EA; it has numerous attractive advantages to solving constrained optimization problems quickly because it is implemented, includes few control parameters, and achieves top rank in many computations 28 .Numerous DE variants have been applied in the literature to find solutions to constrained-type engineering problems.In this work, a constrained composite.DE (C2oDE) global optimizer 25 is proposed and added with two different CHTs to find the balance between constraints and objective functions.The framework of the proposed C2oDE optimization algorithm is introduced in the next section.

C2oDE
In the C2oDE algorithm, differential vectors generate offspring 29 .Fundamentally, there are four stages in the proposed algorithm, in the first stage randomly generation of the initial population − → x t i (i ∈ {1 . . .NP}) in the range of lower and upper bound of search space.After that in the second stage, mutation operators are used for the generation of mutant vector − → v t i (i ∈ {1 . . .NP}) , in this stage three type of mutation operators were used, such as. 1) current-to-rand/l The random differentiation shows the best solution for current generation t and rand.Each mutation vector has distinct features ; for example, the mutation vector given in Eq. ( 30) can explore the entire search space and increase diversity.However, Eqs. ( 31) and (32) accelerate the convergence to get information from the best individual.In the third step trial vector − → u t i is generated using a binomial crossover operator between each pair of − → v t i and − → x t i described as: where, x t i,j , u t i,j and v t i,j are the j th dimension of − → Correspondingly, CR is the rate of crossover, and J rand is the integer number randomly produce between 1 to D. Finally, in the fourth step the selection operator is applied among the x t i and u t i to find the candidate for the next population using Eq.(34) and Fig. 5 shows the framework of the proposed C2oDE algorithm.It can be noticed from Fig. 5 that, for each target vector three off-springs are generated with distinct advantages of exploration and exploitation using trail vector generation strategy and pool of parameters.However, OPF problems are constrained optimization problems; therefore, there must be a compromise between objective function and constraint.Therefore, to balance constraint and objective function, two different CHTs are incorporated in this work at the phase of preselection and selection, as shown in Fig. 5.As stated in No Free Lunch (NFL) 20 , using various CHTs rather than single ones at different stages of EAs is better.Thus, the feasibility rule (FR) and ε constrained method (ECM) two CHTs are implemented with the proposed algorithm at the preselection phase and selection to select feasible trial vectors and populations for the next generation, respectively.OPF problems are very highly complicated.Therefore, a restart technique is used to avoid trapping into local optima, and it is triggered when the standard deviation of both the objective function or constraint violation is less than the assigned threshold value.The flow diagram of the proposed C2oDE-FR-ECM algorithm is given in Fig. 6.C2oDE maintains a population consisting of NP target vectors: , and their degree of constraint violation:

Apply Optimal Load Flow (OPF) to compute objective functions and constraint violation
Feval=Feval + NP

Results and comparison
Various standard IEEE power system test networks were selected to judge the effectiveness of the proposed C2oDE algorithm.These include 30, 57, and 118-bus networks applying two different constraint handling techniques (CHTs) at various stages.Table 7 summarizes the parameters of the proposed algorithm for the simulation of standard IEEE networks provided that values of F and CR are [0.8,1.0, 1.0] and [0.2, 0.1, 0.9], respectively.

Comparison among proposed chts
The C2oDE algorithm is compared and tested with the two most widely used CHTs, FR and ECM, at different places, such as at the preselection stage (to select the best trial vector) and selection (population for the next generation).Table 8 presents the statistical values over the 25 independent runs for the individual events of 1 to 14 using FR and ECM constraint handling methods.The columns of Table 8 show the best, mean, worst, and standard deviation of each event over 25 runs.Table 8 indicates that a single method cannot deliver the best statistical results in all the events.Therefore, this paper includes proposed CHTs in two stages to find a feasible trial vector and population for the next generation.Four different C2oDE variants were implemented considering two CHTs at different locations, i.e.C2oDE-FR, C2oDE-ECM, C2oDE-FR-ECM, and C2oDE-ECM-FR.Specifically, in C2oDE-FR and C2oDE-ECM, only the feasibility rule and ɛ constraint method were utilized for the best trail vector and population of the next iteration.However, in C2oDE-FR-ECM, the feasibility rule was used for finding the best trail vector, and ɛ constraint method was used to select the population for the next iteration while in C2oDE-ECM-FR, ECM for the trial vector, and FR was used to select candidates for the next iteration.
The bold numbers shown in Table 8 are the best objective function values in a particular event obtained by methods.Furthermore, in Table 8, C2oDE-FR and C2oDE-FR-ECM outperform compared to C2oDE-ECM and C2oDE-ECM-FR.In contrast, C2oDE-ECM cannot beat any other variant in any study event, whereas C2oDE-ECM-FR only performs better in event1 and 4. On the other hand, FR and FR-ECM obtain the best fitness value, almost an equal number of events.Hence, selecting the proper CHTs for an OPF problem of various events is challenging because the objective function and constraints of OPF are non-linear.On the other hand, C2oDE-FR-ECM has the benefit of converging with the help of FR and exploring the entire search space to get better diversity with the help of ECM.Thus, the combination of FR and ECM at the different phases of the search algorithm, i.e., in C2oDE-FR-ECM, would attain the best value of the objective function or be close to the best fitness in most of the events.The subsequent subsections analyze and discuss the best results according to the objective functions of all the IEEE test systems.

IEEE 30-bus test system
Table 9 shows the results of 30-bus system decision variables (i.e., state and control variables of event 1 to event 10).Column 2 and 3 of Table 9 displays the operating range of decision variables and in all the events, the results of these variables are within their allowable range and give the best value(s) of fitness considering one of the four proposed algorithms.In this work, the generator's output power in MW at the swing bus (P G1 ) and the MVAr rating of all the generators are considered the control variable and treated as inequality constraints during the optimization process.The allowable range of reactive power for all the generators is taken from MATPOWER 26 .Furthermore, simulation results obtained using the three variants of C2oDE by applying CHTs are presented in Table 10 (for single objective) and Table 11 (for multi-objective) compared with the recent methods of similar studies in the literature.Obtained results of proposed CHTs in which all the decision variables (dependent and independent) and constraints are within desirable limit however, in the approach of static penalty, some of these variables are violated and are highlighted with footnotes as shown in Table 10 and.
Table 11.During the optimization process, voltages at the PQ, buses are often found critical, such as near the upper limit (0.95-1.05 p.u). Frequently, failure of the power system components appears due to overvoltage, and it is highly undesirable.
On the other hand, the Voltage deviation (VD) of the IEEE 30-bus system would be 1.2 p.u (24 × 0.05) if the value of voltage level at all load buses is under the permissible limit.However, in the literature in many cases, VD is more than a permissible specified value and is also highlighted with footnotes, as shown in Table 10 and.11.Furthermore, the main goal of this work is to prove effectiveness by merely considering statistical results and establishing the strict agreement of system constraints using various CHTs.It is noticed from Table 10 that values of the objective function in event 1 using FR and ECM-FR give 800.411290$/h and 800.411384, respectively, satisfying all the inequality constraints.In event 2 C2oDE-FR finds the minimum cost of 646.40111 $/h among various CHTs considering the multi-fuel effect however, in event 3 in which fitness function is considered to minimize the maximum L-index (L max ) of PQ buses, FR-ECM obtained the best simulation result of 0.13628 in comparison to the algorithms of past studies.In event 4, the minimization of emissions in (t/h) is 0.204817, almost the same in all the proposed CHTs.Also, the algorithms reported in the literature include SF-DE 6 , MSA 8 , and ARCBBO 12 whereas, in event 5 minimization of active power losses, FR-ECM and FR give the best results of 3.08391 MW and 3.08392 MW compared with the other techniques shown in Table 10.Voltage waveforms of the 30-bus network are given in Figs.7 and 8 and it shows that the output value of voltage (p.u) is within the range of minimum and maximum value without violating any of the constraints.
The fitness function of fuel cost minimization considering valve-point loading proposed in event 6, in which C2oDE-FR obtained the best result of 832.0700 $/h is high compared to basic fuel cost in event 1.However, in events, 7-10 weighted sum multi-objective optimization of various functions is proposed in which the combined effect of various single objective functions decides the output results of optimization algorithms.For example, in event 7, a higher weight charge of fuel cost was preferred to minimize more fuel cost than power loss.Table 11 shows that the single algorithm FR, FR-ECM, or ECM-FR is not able to find the best value of fitness in all the events.In event 7, C2oDE-FR gives the global minimum of combined fitness of 1040.11188compared to other methods.Furthermore, in events 8 to event 10, obtained values of combined multi-objective functions are minimal in FR-ECM compared to all the algorithms.Furthermore, the convergence curve of C2oDE using two CHTs at different phases for events 1, 2, and 6 considering fuel cost as the objective function are indicated in Figs. 9, 10, 11, respectively.
Among the different CHTs, the convergence speed is not strangely different, though rapid and smooth convergence is observed in both FR and FR-ECM.Figures 12, 13, 14 give the convergence curve of event 3 to event 5, respectively.In event 3, the voltage stability index indicator is scrutinized in the fitness function in which the convergence curve is uneven because of the nature of the objective function.Moreover, Figs. 15, 16, 17, 18 show the convergence curve of multi-objective optimization events.The convergence curve of only the best fitness value of CHTs is shown in Figs. 15, 16, 17, 18 for clear visibility and the irregularities between objective functions during convergence due to non-linear relationships among the fitness and independent variables.

IEEE 57-bus test system
The solution of decision variables (i-e.dependent and control) of the 57-bus network and the simulation results of best objective functions among all the methods are demonstrated in Table 12.
Table 12 clearly shows that the decision variables are within the desirable range.However, Table 13 compares all CHTs (FR, FR-ECM, ECM-FR) with the recent literature methods.Minimum and maximum values of a few generators' MVAr ratings are relatively narrow and taken from 26 even though proposed CHTs dully satisfied the generator reactive power limit.Further, the IEEE 57-bus system consists of 50 PQ buses, and the voltage level of feasible solutions of these buses must be within [0.94 to 1.06] p.u range and the cumulative VD would be (50 × 0.06) 3 p.u.The values of VD found to be more than 3 p.u are marked with a footnote in one reference in which static penalty function is used as CHTs.In events 12 and 14 among four CHTs, the results of C2oDE-FR are best; on the other hand, in events 11 and 13 C2oDE-FR-ECM outperformed among all the proposed CHTs, providing all the constraints are within feasible search space.In most of the events, according to the minimization of objective functions of IEEE 57-bus systems, C2oDE-FR and C2oDE-FR-ECM outperform in comparison to the methods of past studies.In event 11 the best value of the objective function is 41,666.2413$/h, the lowest values comparison to the methods as shown in Table 13 also the power loss (14.86981151MW) is best compared to the method available in the literature.Event 12 is the multi-objective, considering fuel cost and VD by C2oDE-FR is 41,774.422,which is close to the value given by SP-DE 6 .Further, multi-objective voltage stability and fuel cost are considered in event 13, in which fuel cost (41,694.089)seems better compared to event 12.In event 14 C20DE-FR outperformed APFPA 13 , however with the expense of fuel cost compared to the previous study.Larger values of shunt VAR compensators (30 MVAr) have been seen in 16 and 7 , hence comparison with the present study is not valid.Figure 19 shows the voltage profiles of the best solution among the different CHTs of event 11 to event 14 for the 57-bus system.
Figure 19 clearly shows that the operating value of the voltage at all the buses is within the minimum and maximum range, such as satisfying voltage constraints so that no bus experiences overvoltage, whereas, in some buses, the voltage level is close to the upper bound.Figure 20 shows the convergence curves of applied CHTs.As  We can notice from the above figures that the convergence curves of voltage deviation and L-index are nonsmooth in multi-objective events.Figure 23 shows the convergence curve of single objective optimization (voltage deviation) in which all methods need many fitness function evaluations to seek the global optimum solution because of the non-linear relation between bus voltage and independent variables in the 57-bus tests network.

IEEE 118-bus test system
Generally, for an increased number of variables, the performance of C2oDE-FR-ECM is found to be superior.Hence, in the large-scale 118-bus test network, an effective combined FR-ECM constraint technique is proposed to show the superiority and scalability of the proposed algorithm.Furthermore, the minimization of basic fuel     cost (Event 15) and real power losses (Event 16) is considered the objective functions of this system.Table 14 shows the calculated parameters and control variables of the best solution found using C2oDE-FR-ECM.
Allowable values of MW and MVAr rating of generators, the voltage level of transformers, and the MVAr rating of shunt VAR compensators are taken from 26 , and Table 14 clearly shows that in events 15 and 16, all the control variables are fully satisfied the minimum and maximum limit.The results of event 15 and event 16 are that the basic fuel cost is 134,943.8$/h and active power losses are 16.79906MW, respectively.Figure 24 shows the voltage profile of all the buses and the minimum and maximum limits, while Fig. 25 gives the convergence curves of events 15 and 16.
However, Table 15 shows the comparative results of the proposed C2oDE-FR-ECM with the recently implemented DE variants in the literature.Table 15 shows that the proposed algorithm finds a better approximate optimal solution than all the other state-of-the-art evolutionary algorithms.

Conclusion
Optimal power flow (OPF) is a highly complex, constrained, and non-linear problem in a power system.In the solution of OPF problems without using suitable CHTs, the decision variable of the system may be violated and given poor safety, ill-functioning protective devices, and unnecessary power losses, especially with a static penalty function.Therefore, during the operation of the power system, constraint handling techniques (CHT) are responsible for optimizing objective functions subject to decision variables, and constraint functions should be within safe limits.Therefore, the application and usefulness of two CHTs, such as feasibility rule (FR) and ε constraint method (ECM), and their combinations with outstanding global optimizer C2oDE (C2oDE-FR, C2oDE-ECM, C2oDE-FR-ECM, C2oDE-ECM-FR) have been presented and used to solve OPF problem taking into various non-linear constraints.Three standard test networks, small to large-scale power system networks such as IEEE 30, 57, and 118-bus, are scrutinized to solve OPF problems with the CHTs group that helps achieve the best feasible solution in most of the events.A comparative analysis of the four techniques reveals the challenge    nearly optimal solutions in most events.However, it does not guarantee the most optimal solution or rapid convergence in all events.Nonetheless, the significance of an efficient constraint-handling technique cannot be overstated.As our study demonstrates, inadequate CHT, mainly the penalty approach, may unknowingly lead to violations of network parameter limits.Therefore, to ensure a feasible solution to the OPF problem, the power system constraints must be within defined limits, is essential for its secure and proper functioning.The recommended configurations of CHT effectively bring the network to the desired state compared with several other methods in the past study.

Figure 3 .
Figure 3. Output power vs fuel cost of single and multi-fuels.

2 )
Modified rand-to-best/l 3) current-to-best/l where, − → x t r1 to − → x t r5 are the mutually different decision vectors randomly selected from 1 to NP individuals, − → x t best

Figure 7 .
Figure 7. Event-1 to event-5: Voltage profile of best solution For IEEE 30-bus systems.

Figure 8 .
Figure 8. Event-6 to event-10: Voltage profile of best solution For IEEE 30-bus systems.

Table 1 .
Summary of IEEE 30-bus test system under study.

Table 2 .
Coefficients of cost and emission (generators) for 30-bus network.

Table 3 .
Multi-fuel cost coefficients of generators 1 and 2 of the IEEE 30-bus test system.

Table 4 .
Data of IEEE 57-bus network under study.

Table 7 .
PARAMETERS OF Proposed C2oDE ALGORITHM.

Table 8 .
Statistical summary of FR and ECM of event 1 to event 14.

Table 9 .
Simulation results of event 1 to event 10 considering the best algorithm for a 30-bus network.

Table 10 .
Comparison of results of proposed algorithms with the past studies of the 30-bus single objective.
a Voltage level at the PQ bus is violated.

Table 11 .
Comparison of results with the recent methods of proposed algorithms for 30-bus multi-objective.

Table 12 .
Simulation results of the best algorithm for 57-bus network.

Table 13 .
Comparison of proposed algorithms with the past studies of IEEE 57-bus system.a Voltage on the PQ bus is violated, the solution is infeasible.b Considered the large limits of shunt VAR compensators.

Table 15 .
Comparison of proposed algorithms with the past studies of IEEE 118-bus system.